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Abstract We prove the continuity and we give a Holder estimate for the value function as- 
(T) | sociated with the L 1 cost of the control-affine system q = fo(q) + J2T=i u jfj(l)^ satisfying the 

strong Hormander condition. This is done by proving a result in the same spirit as the Ball-Box 
theorem for driftless (or sub-Riemannian) systems. The techniques used are based on a reduction 
of the control-afhne system to a linear but time-dependent one, for which we are able to define a 
generalization of the nilpotent approximation. Finally, we also prove the continuity of the value 



function associated with the L 1 cost of time-dependent systems of the form q = Y^jLi u jfj(l)- 
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1. Introduction 
A sub-Riemannian control system on a smooth manifold M is a control system in the form 



u 

o 

^^ in 

(1) 7(*)=£«i(*)/i(7(*)), a.e.t€[0,T], 

where u : [0,T] — )■ R m is an integrable control function and {/i, • • • , /m} is a family of smooth 

vector fields satisfying the Hormander condition, i.e., such that its iterated Lie brackets generate 

the whole tangent space at any point. The length of a curve 7 solving (TXJ) , is then defined as the 

_j_ . length(7) = min J Q \u\ dt, where the minimum is taken over all the possible u(-) satisfying the 

\Q ' above ODE. Due to the linearity of the system w.r.t. u, this length will be independent of the 

VO . parametrization of 7. Finally, we define 

dSB.(q, q') = inf {length(7) : 7 : [0, 1] -> M is a solution of (TTJ), 7(0) = q and 7(T) = q'}. 

By the Hormander condition, dsn is a distance, called Carnot-Caratheodory distance, endowing 
M with a natural metric space structure. A manifold considered together with a sub-Riemannian 
control system is called a sub-Riemannian manifold. 

Define A 1 = span{/i, . . . , f m } and A s+1 = A s + [A*, A 1 ], for every s £ N. Under the 
hypothesis that {/i,...,/ m } is equiregular, i.e., that, for each s G N, the dimension of A s (q) 
is independent of q G M, the Hormander condition implies that there exists a (minimal) r G N 
such that A r (q) = T q M for all q G M. Such r is called the degree of non-holonomy of the 
sub-Riemannian control system. 
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A fundamental result in the theory of sub-Riemannian manifolds is the celebrated Ball-Box 
theorem (see for example (5J). This theorem gives a rough description of the infinitesimal shape 
of the sub-Riemannian balls. Namely, at any point q of an equiregular sub-Riemannian manifold, 
the sub-Riemannian ball of small radius e is equivalent, in privileged coordinates, to the box 

[— e, e] x ... x [— e, e] x ... x [-e s ,e s ] x . . . x [-e s , e s ] x . . . x [-e r , e r ] x . . . x [~s r , e r ] . 

V v ' V v ' V v ' 

dimA 1 dim A s — dim A s_1 dim A r — dim A 7 " -1 

By this we mean that, for some constant C > 0, the sub-Riemannian ball is contained in a box 
of side Ce, and contains a box of side e/C. 

This fact has a plethora of applications. First, it allows to prove a Holder regularity estimate 
with respect to the Euclidean distance in coordinates, namely that, for q' sufficiently close to q, 
it holds 

(2) \q'\<d SR (q,q')<\q'\^. 

Here we use "<" to denote an inequality up to a multiplicative constant, independent of q' . Then, 
among many others, it is a fundamental step in the computation of the Hausdorff dimension of 
the manifold (see [21]), and it is used to obtain asymptotic estimates on the heat kernel (see e.g., 
[231 1201 H3J [2]). Moreover, it is the main tool in computing the asymptotic equivalents of the 
entropy and the complexity of curves (see e.g., [HI HH [22j [T5l H6] ). 

In this paper, we focus on a very important generalization of the control system ([T]), namely 
on control- affine systems. These systems are obtained by adding to ([1]) an uncontrolled vector 
field /o, called the drift, and are in the form 

m 

(3) 7(t) = / (7(*)) + J2 u ^ /<(7(*)). a -e- t £ [0, T]. 

These kind of systems appears in plenty of applications. As an example we cite, mechanical 
systems with controls on the acceleration (see e.g., [TUJ, [Z]), where the drift is the velocity, or 
quantum control (see e.g., [T2], [5]), where the drift is the free Hamiltonian. We always assume the 
strong Hormander condition, i.e., that the family {/x, . . . , f m } satisfies the Hormander condition. 
The cost of a curve 7 solving ([3]) is Cf a {l) — mm Jo \ u \ dt. Unlike the sub-Riemannian length, 
due to the presence of the drift, the cost depends on the parametrization of the curve. Finally, 
the value function, between q, q' £ M, of the control system at time T > 0, is defined as 

p f T °(q,q') = inf{c /o ( 7 ): 7 : [0,T'] -> M solves ©, 7(0) = q, j(T') - q', and T' < T}. 

Assume now that the drift is regular, in the sense that there exists s £ N such that fo(q) £ 
A s (g) \ A s ^ 1 (q), for any q £ M, where A s is defined through the vector fields {/1, . . . , /,„} as 
before. Our main result is, then, a generalization of (J5J to this context. 

Theorem 1. Let z = (z\, . . . , z n ) be a system of privileged coordinates at q for {/1, . . • , fm} , 
rectifying /o as the k-th coordinate vector field d Zk , for some 1 < k < n. Then, for sufficiently 
small T and e, if Pj!(q, q') < e then in z coordinates it holds 

dist (z(q'),z(e [0 ' T1/o g)) < p f 4\q,q') < dist (z(g')^(e [0 ' T1/o (7)) ^ . 

Here for any x £ M. n and A C W 1 , dist(a;, A) — mi ye A \x — y\ denotes the Euclidean distance 
between them, e*^° denotes the flow of fo, and r is the degree of non-holonomy of the sub- 
Riemannian control system defined by {/1, . . . , /,„}. 

In this result, instead of the Euclidean distance from the origin that appeared in ([2]), we 
have the distance from the integral curve of the drift. This is due to the fact that moving in 
this direction has null cost. As in the sub-Riemannian case, Theorem Q] is a consequence of 
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an estimate on the shape of the reachable sets, contained in Theorem [2JJ Moreover, although 
Theorem Q] seems a natural generalization of ©, the shape of the reachable sets described in 
Theorem [57] is much more complicated than the boxes of the sub-Riemannian case, yielding a 
more difficult proof. Theorem [1] and [27] represent the key step for generalizing the estimates on 
the complexity of curves from sub-Riemannian control systems to control- affinc systems. 

It is worth to mention that these results regarding control-affine systems are obtained by 
reducing them, as in [3J, to time-dependent control systems in the form 

m 

(4) 7(<) =$>*(<) #(7(*)), a.e. ie[0,T], 

4=1 

where /* = (e~ f ^°)*/i is the pull-back of fa through the flow of the drift. On these systems, that 
are linear in the control, we are able to define a good notion of approximation of the control vector 
fields. Namely, in Section [33] we will define a generalization of the nilpotent approximation, used 
in the sub-Riemannian context, taking into account the fact that in the system ([¥]), exploiting the 
time, we can generate the direction of the brackets between /o and the fjS. This approximation 
and an iterated integral method yield the correct estimates on the reachable set, contained in 
Theorem [H 

The paper is divided in three sections. In Section [2] we recall some generalities and definitions 
regarding sub-Riemannian control systems, used in the following sections. In Section [3J we 
consider control systems in the form (jj), and we prove the continuity of the value function 
for general time-dependent vector fields. Then, in Theorem 1191 restricting then to the case 
where the time dependency is explicitly given as ft = (e _ *^°)*/i, we establish some estimates 
on the reachable sets, in the same spirit as the Ball-Box theorem. Finally, in Section [4] we 
consider control-affine systems. After proving the relation between control-affine systems and 
time-dependent systems, we prove the continuity of the value function. Then, in Lemma 1281 
exploiting the affine nature of the control system, we give an upper bound on the time needed to 
join two points q and q' as a function of pl^ (q, q 1 ). From this fact and the estimates of Section [3] 
Theorems l2"ol and |2"T1 follow. Theorem Q] is then a particular case of Theorem [2^1 that holds under 
slightly milder assumptions on /o and {/i, . . . , f m }. 

2. SUB-RlEMANNIAN GEOMETRY 

Throughout this paper, M is an n-dimensional connected smooth manifold. In this section 
we recall some classical notions and results of sub-Riemannian geometry. 

2.1. Sub-Riemannian control systems. A sub-Riemannian (or non-holonomic) control sys- 
tem on M is a control system in the form 

III 

(SR) i = ^Uifi{q), q€M, u = (u x , . . . ,u m ) € R m , 

where {/i, . . . , f m } is a family of smooth vector fields on M. We let f u = YlTLi u i fi- 

An absolutely continuous curve 7 : [0,T] —> M is (|SR[) -admissible if there exists a control 
u e L 1 ([0,T],R m ) such that j(t) = f u (t)(l(t)), for a.e. t € [0,T]. The curve is said to be 
associated to any such control. The length of 7 is defined as 

(5) length(7) = min |M| l i([o,t],r™), 

where the minimum is taken over all controls u such that 7 is associated with u. It is at- 
tained, due to convexity reasons. Notice that, by definition, length^) is invariant under time 
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reparametrization of the curve. The distance induced by the sub-Riemannian system on M is 
then defined as 

dsn(q, q) = inf {length^) : 7 (jSRP -admissible and 7 : q ~~» </}, 

where 7 : q -^ q' stands for 7 : [0, T] — > M, for some T > 0, 7(0) = q and 7(T) = q' ■ 

Let A be the C°°-module generated by the vector fields {/i,...,/ m } (in particular, it is 
closed under multiplication by C°°(M) functions and summation). Let A 1 = A, and define 
recursively A s+1 = A s + [A S ,A], for every s G N. Due to the Jacobi identity A s is the C°°- 
modulc of linear combinations of all commutators of /1, . . . , f m with length < s. For q G M , let 
A s (q) = {f(q) '■ f G A s } C T q M. We say that {/1, . . . , / m } satisfies the Hormander condition 
(or that it is a bracket-generating family of vector fields) if U s >i A s (g) = T q M for any q G M. 
In the following we will always assume this condition to be satisfied. 

By the Chow-Rashevsky theorem (see for instance [T]), the hypothesis of connectedness of 
M and the Hormander condition guarantee the finiteness and continuity of dsR with respect to 
the topology of M. Hence, the function dsR, called sub-Riemannian or Carnot-Caratheodory 
distance, induces on M a metric space structure. The open balls of radius e > and centered at 
q G M, with respect to ^sr, are denoted by Bsr(<7, e)- 

We say that a fSRj-admissible curve 7 is a minimizer of the sub-Riemannian distance between 
q, q' G M if 7 : q -^> q' and length^) = dgYi(q,q')- Equivalently, 7 is a minimizer between 
q,q' G M if it is a solution of the free-time optimal control problem, associated with (jSRI) . 



(6) I|m||li(o,t) , , 

/o \ 



J2tf(t)dt ^min, 7(0) = q, j(T)=c{, T > 0. 

3=1 



Indeed, the sub-Riemannian distance is the value function associated with this problem. It is a 
classical result that, for any couple of points q, q' G M sufficiently close, there exists at least one 
minimizer. 

Remark 2. The optimal control problem ^ is equivalent to the following, with p > 1 and 
T > fixed, 




(7) HU»(o,T)= / W\ p dt) -Main, 7 (0) = q, j(T) = q', 



In fact, due to the invariance under time reparametrization of system (|SR|) . in © we can fix 
cither T > or the Euclidean norm of u. Moreover, by the Cauchy-Schwarz inequality, for any 
p > 1, letting p' be the conjugated exponent to p (i.e., 1/p + 1/p' = 1), we get IHIl^o.t) < 
r 1 / p '||-u|| i p(o,T) J with the equality holding if and only if |u| is constant. From these two facts, it 
is easy to check that minimizers of the optimal control problem (0) coincide with the minimizers 
of ([B]) with constant norm. Thus the two optimal control problems are equivalent. 

Remark 3. This control theoretical setting can be stated in purely geometric terms. Indeed, it 
is equivalent to a generalized sub-Riemannian structure. Such a structure is defined by a rank- 
varying smooth distribution and a Riemannian metric on it (see pQ for a precise definition). In 
a sub-Riemannian control system, in fact, the map q n* span{/i(g), . . . , f m (q)} C T q M defines a 
rank-varying smooth distribution, which is naturally endowed with the Riemannian norm defined, 
for v G A(q), by 



(q, v) = inf I \u\ = sju{ + h u 2 m : f u (q) = v 
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The pair (A, g) is thus a generalized sub-Riemaimian structure on M. Conversely, every rank- 
varying distribution is finitely generated, see [51 [Tl HI 113). and thus a sub-Riemannian distance 
can be written, globally, as the value function of a control system of the type (jSRI) . 

2.2. Privileged coordinates and nilpotent approximation. We now introduce the equiva- 
lent, in the sub-Riemannian context, of the linearization of a vector field. This classical procedure, 
called nilpotent approximation, is possible only in a carefully chosen set of coordinates, called 
privileged coordinates. 

Since {/i, . . . , f m } is bracket-generating, the values of the sets A s at q form a flag of subspaces 
of T q M, 

A 1 (q)cA 2 (q)c...c£ r (q)=T q M. 
The integer r — r(q), which is the minimum number of brackets required to recover the whole 
T q M, is called degree of non-holonomy (or step) of the family {/i,...,/ m } at q. Set n s {q) — 
dim A s (q). The integer list (ni(q), . . . ,n r (q)) is called the growth vector at q. From now on we 
fix q G M, and denote by r and (rix, . . . , n r ) its degree of non-holonomy and its growth vector, 
respectively. Finally, let w% < . . . < w n be the weights associated with the flag, defined by Wi = s 
if n s _i < i < n s , setting uq = 0. 

For any smooth vector field /, we denote its action, as a derivation on smooth functions, by 
/ : a G C°°(M) n- fa G C°°(M). For any smooth function a and every vector field / with / ^ 
near q, their (non-holonomic) order at q is 

ord q (a) = min{s G N: 3i u ...,i s G {l,...,m} s.t. (f h ...f is a)(q) ^ 0}, 

ord q (/) = maxjcr G Z: ord g (/a) > a + ord g (a) for any a G C°°(M)}. 

In particular it can be proved that ord g (a) > s if and only if a(q') = 0(dsYi(q' ,q)) s ■ 

Definition 4. A system of privileged coordinates at q for {/i,...,/m} is a system of local 
coordinates z — (z\, . . . , z n ) centered at q and such that ord g (zi) = Wi, 1 < i < n. 

For any point q G M there always exists a system of privileged coordinates around q. Consider 
such a system z = (z\, . . . , z n ). We now show that this allows to compute the order of functions 
or vector fields in a purely algebraic way. Given a multiindex a — (ai,...,a„) we define 
the weighted degree of the monomial z a = z" 1 • • • z"™ as w(a) = w\ai + • • • + w n a n and the 
weighted degree of the monomial vector field z a d z . as w(a) — Wj . Then one can prove that, given 
a G C°°(M) and a smooth vector field /, with Taylor expansion 

a ( z ) ~ X!«aZ a and f(z) - ^f a jz a d Zj , 

their orders at q can be computed as 

ord g (a) = min{u>(a) : a a ^ 0} and ord q ( f) = m.m{w(a) — Wj : fa.jy^O}. 

A function or a vector field are said to be homogeneous if all the nonzero terms of the expansion 
have the same weighted degree. 

We recall that, for any a, b G C°°(M) and any smooth vector fields f,g, the order satisfies the 
following properties 

ord g (a + b) = min{ord g (a), ord g (&)}, ord g (a6) = ord g (a) +ord q (6), 

ord q (f + g) =min{ord g (/),ord g ( 3 )}, ord g ([/,.g]) > ord q (/) + ord g (5). 

Consider fi, 1 < i < m. By the definition of order, it follows that ord ? (/i) > —1. Then we 
can express fi in coordinates as 

z*fi = ^2 (hij +rij)d Zj , 



HOLDER CONTINUITY OF THE VALUE FUNCTION FOR CONTROL-AFFINE SYSTEMS 6 

where z* is the push-forward operator on vector fields associated with the coordinates, defined 
asz,/ = dzo/oz _1 , hij are homogeneous polynomials of weighted degree Wj — 1, and r^ are 
functions of order larger than or equal to Wj . 

Definition 5. The nilpotent approximation at q of fi, 1 < i < m, associated with the privileged 
coordinates z is the vector field with coordinate representation 

n 
3=1 

The nilpotentized sub-Riemannian control system is then defined as 

m 

(NSR) ? = £>(*)£• (g). 

The family of vector fields {/i, . . . , f m } is bracket-generating and nilpotent of step r (i.e., every 
iterated bracket [fi 1 ,[..., [fi k _ 1 , fi k })} of length larger than r is zero). 

The main consequence of the nilpotent approximation is the following (see for example [5J 
Proposition 7.29]). 

Proposition 6. Let z = (zi,...,z„) be a system of privileged coordinates at q G M for 
{fi,...,f m }. For T > and u e ^([0, T]; R m ), with \u\ = I, let 7 (-) and *{{■) be the tra- 
jectories associated with u in (|SR[) and (INSR[) . respectively, and such that 7(0) = 7(0) = q. 
Then, there exist C, To > 0, independent of u, such that, for any t < Tq, it holds 

(9) \z l ( 1 (t))-z i m))\<Ct w >+\ i=l,...,n. 

We recall, finally, the celebrated Ball-Box Theorem, that gives a rough description of the 
shape of small sub-Riemannian balls. 

Theorem 7 (Ball-Box Theorem). Let z — (zi, . . . ,z n ) be a system of privileged coordinates at 
q E M for {/1, . . . , f m }- Then there exist C,eo > such that for any e < Eq, it holds 

Box ( —e j C B SR (q, e) C Box (Cs) , 

where, -Bsr((J, e) is identified with its coordinate representation z(Bsn(q 1 e)) and, for any r\ > 0, 
we let 

(10) Box(ri) = {zem. n : \zi\<r) Wi }, 

Observe that the first inclusion follows directly from the definition of privileged coordinates. 
As a corollary of the Ball-Box Theorem, we get the following result on the regularity of the 
distance. 

Corollary 8. Let z = [z\, ■ ■ ■ , z n ) be a system of privileged coordinates at q £ M for {/1, . . . , f m }. 
Then there exists C, e > such that 

^W)\ < d SR (q,q') < C\z(q')\ 1/r , q' € B SR (q,e). 

3. Time-dependent systems 

3.1. Time-dependent control systems. We now consider a more general situation. Namely, 
we consider on M the time-dependent non-holonomic control system 

m 

(TD) q = J^Uift(q), q € M, u = ( Ul , . . . , u m ) e R m , tel, 
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where / = [0, b) for some b < +00 and {/', . . . , /^} is a family of non-autonomous smooth vector 
fields, with smooth dependence on the time parameter. We let /* = Yl^Li u i ft- 

In analogy with the autonomous case, we define (ITD[) -admissible curves as absolutely contin- 
uous curves 7 : [0,T] C / — 5- M such that j(t) = f^r t )(j(t)) for a.e. t £ [0, T], for some control 
u £ L 1 ([0,T],1R" 7 '). Observe, however, that contrary to what happens in the sub-Riemannian 
case, the (jTDp -admissibilitv property is not invariant under time reparametrization, e.g., a time 
reversal. Thus, we define the cost (and not the length) of 7 to be 

c(j) =min||u|| L i([o > T],R"») ) 

where the minimum is taken over all controls u such that 7 is associated with u and is attained 
due to convexity. The value function induced by the time-dependent system is then defined as 

p(q, q') = inf {0(7) : 7 is (|TDI) -admissible and 7 : q ~-» q'}. 

Clearly, the value function is non-negative. It is not a metric since, in general, it fails both to be 
symmetric and to satisfy the triangular inequality. Moreover, as the following example shows, p 
could be degenerate. Namely, it could happen that q ^ q' but p{q 1 q') = 0. 

Example 9. Let M = R, with coordinate x and consider the vector field /* = (1 — t)~ 2 d x 
defined on [0, 1). For any xo £ 1, in ^ 0, and for any sequence t n f 1, let u n £ L 1 ([0,i„]) be 
defined as u n = (1 — t n )xo. By definition, each u n steers the system from to xq. Hence, 



Pi(0,x o ) < inf ||w„||Li([o,t„l) = inf / (I - t n )x dt = x in£t n (l-t n ) =0. 

nSN "GNJq n£N 

This proves that, for any xq £ R, pi(0 7 xo) = 0. 

For T > 0, q € M and e > 0, we denote the reachable set from q with cost less than e by 

ll(q,E) = {q' eM: p(q,q') < e}. 

We will also consider the reachable set from q in time less than T > and cost less than e, and 
denote it by TZr(q 7 e). Clearly 7£x(<7, e) C 1Z(q,e). 

In general, the existence of minimizers for the optimal control problem associated with (jTDp 
is not guaranteed. We conclude this section with an example of this fact. 

Example 10. Let M = R, with coordinate x, and consider the vector field /* = e~ t d x for 

t e [0,1). Fix xq £ R, x ^ 0. Observe that, for any T > and any control u £ L 1 ([0,T]) 
steering the system from to x$ , it holds 



(11) |so| = 



T 

-t 



tl 



(t)e~*d< 



T 

t 



< / |«(t)|e { dt< ||u|| L i ([ o,T])- 



(i 



This implies p(0,xo) > \xq\. Let now u n £ L 1 ([0, 1/n]) be defined as u n (t) = x^ne 1 . Clearly u n 
steers the system from to xq. Moreover, 

e»-l 

p(0,x o ) < inf ||M„|| L i ([01/rl]) = |x | inf — j — = |x |. 

n 

This proves that p(0,xo) = |xo|. Hence, the non-existence of minimizers follows from (|lip. 

3.2. Finiteness and continuity of the value function. In this section, we extend the Chow- 
Rashcvsky Theorem to time-dependent non-holonomic systems, under the strong Hormander 
condition, whose definition follows. 

Definition 11. We say that a family of time-dependent vector fields {/', . . . , f^tei satisfies the 
strong Hormander condition if {f x , . . . , /^J} satisfies the Hormander condition for any to £ !• 
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As we will see later on in Section 01 when considering families of time-dependent vector fields 
of the form /' = (e - *^°)*/i this condition is equivalent to the strong Hormander condition for 
the affine control system with drift /o and control vector fields {/i, . . . , / m }- 

From now on we will assume that the following holds. 

The family of smooth vector fields {/*,.-., /m}*e I ■> depends smoothly on t 
and satisfies the strong Hormander condition. 
This section will be devoted to the proof of the following. 

Theorem 12. Assume that {/*, . . . ,/4,}te/ satisfies ([Ho). Then, the function p : M x M —$■ 
[0,+oo) is continuous. Moreover, for any to G / and any q,q' G M, letting e?sR be the sub- 
Riemannian distance induced by {fi°, ■ ■ ■ , fm}, it holds p(q, q') < dsR^q, q')- 

Now, we need to introduce some notation. Following [5], the flows between times s, t G R of 
an autonomous vector field / and of a non-autonomous vector field r H> f T will be denoted by, 
respectively, 



,(*-«)/ . M -^ M and e^> f f T dr:M^ M. 



Fix q G M and assume, for the moment, that to = 0. Let £ G N and T — (ii,...,ii) G 
{1, . . . , m} . For any T G /, 7~ > 0, we define the switching end-point map at time T and 
associated with J 7 to be the function E-r.F '■ K f — > M defined as 



(12) 



ErAO = "^ /_, J&fldTo.-.o^J 1 ifc /u dr (g) 
= $$J e _JtfT t r dTO...o£$l a /^ dr (g). 



Here we applied a standard change of variables formula for non- autonomous flows. Let then 

[ixfij ao<r<i/i, 



(13) a T T ,T 



&ft mT ifi/e<r<2/i, 



Ueft (i ' m)T if(/-l)//<r<l, 



so that we can write 

EtAZ) = u$ f 9 T T ,A0dr(q). 
Jo 

Clearly, t \-> ex^/ gly jr(t;) dr (q), t G [0,1], is a (ITDp -admissible trajectory. Thus, E-r,r{£), 

T > 0, is the end-point of a piecewise smooth (|TD|) -admissible curve. 

We recall that, by the series expansion of exp (see [5]), for any non- autonomous smooth vector 

field / T , it holds exp7 * f T dr (q) = e* f ° (q) + 0{t 2 ). Thus, we can define 



f 1 

EoAt) = ^ E T,A0 = & f " °-°e 6 fl (?) = $$ / <?o>(0 dr 
'4-0 Jo 



(«), 



where, 30,^(0 * s defined in (fT5)) . Then 1 1-4 exp/ go^(0 dr (g), t G [0, 1], is an (|SR|) -admissible 
curve for the sub-Riemannian structure defined by {/j 3 , . . . , /^} and £a,.f(0 is the end-point of 
a piecewise smooth trajectory in (ISRI) . 

After [24], we say that a point q' G M is (|TDjl -reac/ta6fe from g at time to = 0, if there exist 
t G N, T G {1, .. . ,to}^, T > and £ G M f , such that E T ^{£) = q'. In this case it is clear 
that p(q,q') < X^i l&l- Moreover, if £' i->- Ef t jr(^') nas rank n at £, the point q' is said to be 
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(|TDp -normally reachable at time to = 0. Finally, the point q' is said to be (ISRI) -reachable or 
(|SR| -normally reachable for the vector fields {/° , . . . ,/^}, if these properties holds for T = 0. 

In the case to > 0, taking T > such that 7~+ to G 7 and changing the interval of integration 
in (fT2|) from [0,7] to [io>^o + 7], it is clear how to define (jTDI) -reachable and (|TD|) -normally 
reachable points from q at time to, and (|SR| -reachable and (JSR)) -normally reachable points for 
the vector fields {/*°, . • • , /,'?}■ 

The proof of the following lemma is an adaptation of [231 Lemma 3.1]. 

Lemma 13. Let q' G M be (ISRI) -normally reachable for the vector fields {f^, ■ ■ ■ , /„} from q, 
by some £ G N, £ G R and J 7 6 {1, . . . , m}^ . Then, there exist £o,To > such that, for any 
e < £o, the point q' is (|TPp -normally reachable at time to, by the same £ and J- , and some 
£' G R e , with £\ |& - ^-| < e, and any T < To- 
Proof. Without loss of generality, we assume to = 0. 

Let [/ C l' be a neighborhood of £ such that Eo^r has still rank n when restricted to it. 
Then, there exists B — {x: £?. |xj — £y| < e} C [7 such that Eq,f maps diffeomorphically 
a neighborhood of B in U onto a neighborhood of q. It follows, from standard properties of 
differential equations, that, for T > sufficiently small, the map E-y-^jr is well defined on B and 
that Eq- : jr — >• -E'c.f as T 4- in the C 1 -topology over B. Thus, there exists 71 > such that, 
for 7" < 71, Eq-^jr has rank n at every point of B. 

Since the map Eo.r is an homeomorphism from B onto a neighborhood of q, and E-j-.jr — > 
Eo.jr uniformly as T I 0, it follows that there exists a fixed neighborhood V of q and T2 > such 
that V C Eq-^jr(B), for any T < T<x- Then, for any T < min{7!,72}, there exists £' € -B such 
that the point q' = -Et..f(£') is (|TD|) -normallv reachable. □ 

We will use the following consequence of Lemma [T31 We remark that the result holds even if 
{/', . . . , /4,}te/ satisfies the Hormander condition only at time to G I. 

Lemma 14. Let dsR fee t/ie sub-Riemannian distance induced by {fi°, ■ . ■ ,fm}> then for any 
t\ G /, such that t\ — to > is sufficiently small, and for any q, q' G M it holds 

inf {0(7) : 7 : [to,ti] — > M is (ITDI) -admissible, j(tp) — q and j(t\) = q'} < d$Yi(q,q')- 

In particular, p(q,q') < dsR(q,q')- 

Proof. Fix e > 0. By Chow's theorem it is clear that q' is fSR|-reacheable from q. Moreover, since 
there exist (|SRp -normallv reachable points from q' arbitrarily close to q 1 (see e.g., [TJ Lemma 
3.21]), follows that q' is always (ISRP -normallv reacheable from q by £ such that ^ |£j| < 
dsn(q, q') + s/2. Hence, by Lemma IT3l if e and r\ > are sufficiently small, we have that q' is 
(|TDp -normallv reachable from q at time to by £' such that J2 ■ |£'| < dsn(q,q') + £ and T < t\. 
This clearly implies that 

inf{c(7) : 7 is (jTDp -admissible. 'v(tn) = q and 7(^1) = q'} < dsR(q, q') + £■ 

Finally, the lemma follows letting e \. 0. □ 

We now prove the main theorem of the section. 

Proof of Theorem \1'SX By Lemma 1141 we only need to prove the continuity of p. We will prove 
only the lower semicontinuity, sinche the upper semicontinuity follows by similar arguments. 

We start by proving the lower semicontinuity of p{q, ■) at q' . Consider a sequence q^ — > q' 
and let u^ G L 1 ([0,Tfc],R m ) be controls such that each one steers system (|TPp from q to qk and 
liminf n p*°(q, qk) = liminf n ||itfc|| L i. Then, by Lemma 1141 for any e > there exists a sequence 
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of Tfc > and a sequence of controls vj- G L 1 ([T/ £ ,Tfc],R m ) all steering system (ITD|) from qk to 
q' and such that Hufcllnfm. f .1 K m ) — ^Sr(9*:, g') + £• Since dsR(gfc, g') - ► 0, this implies that 

P(q,q') < n ^n o f ||«fc||Li([o,T fc ]^m) + IMI L i([T fc ,f fc ],R™])J = liminf />(g,g fc ) +e. 

Letting e 1 proves that p(g, •) is lower semicontinuous at q' . 

In order to prove the lower semicontinuity of p(-,q') at g, let us define 

ip s (p) = inf{c(7) : 7 : [e, T] C / -> M is (|TD[) -admissible and 7 : p — * q'}. 

We claim that for any p G M it holds that y e (p) — >• p(p,q') as e 1 0. Since it is clear that 
fe(-) > p(-, <l'), it suffices to prove that 

(14) lim(p £ (p) < p(p,q') for any p G M. 

To this aim, fix p G M and 77 > and let 7 : [0, T] — > M be such that 0(7) < p(p, q') + rj. It is clear 
that 7(2e) — > p as s 1 0, and hence that /o(p, 7(2e)) -» as s i 0, by the first part of the proof. 
Thus, for any e > sufficiently small, there exists a (ITDp -admissible curve 7 e : [e, 2e] — > M such 
that 7 e : p ~-» 7(2e) and c(j e ) < p(p,j(2e)) + r). By concatenating 7 e with 7|[2e,T]> we get that 

V«(p) < c(7„) + c( 7 ) < p(p, 7 (2e)) + p(p, g') + 2t?. 

Letting e i and then 77 ! 0, this proves (fT4")l and thus the claim. 

Let now qk —> q and fix rj > 0. By Lemma 1 141 this implies that p(qk,q) — > and that for 
any e > sufficiently small, there exists a ((TDj-admissible curve 7 e : [0, s] — > M such that 
7e : g fc ~> g and c(7 e ) < p(qk,q) + V- Hence 

p(qk,q') < c(7 e ) + y e (g) < p(g fc , g) + y? 6 (g) + »?• 

By the previous claim, letting £,77 J, 0, this implies that p(qk,q') < p(qk,q) + p{q,q')- Since 
p(qk,q) —^ 0, taking the liminf as fc — > +00, this proves the lower semicontinuity of p(-,q') at g, 
completing the proof. □ 

Remark 15. From the proof of Theorem [T^l it follows that hypothesis QHqD is not sharp. Indeed, 
the following is sufficient to prove the theorem. 

The family of smooth vector fields {/', . . . , /^} t6 /, depends smoothly on t, 
(Hi) and satisfies the strong Hormander condition at t — and in an open 

neighborhood of sup /. 

We will conclude this section by showing that, in our framework, it is necessary to assume the 
Hormander condition on both ends of I. Although outside the scope of the present work, we 
remark that stronger assumptions on the regularity of the vector fields, i.e., that they are uni- 
formly Lipschitz, would allow to prove Theorem [T21 assuming only that {/*, . . . , f^}tei satisfies 
the Hormander condition at one time to G /. 

The following example proves that if the family {/', . . . , ffnjtei satisfies the Hormander con- 
dition only near t = 0, then the value function is in general not continuous. Through a slight 
modification, the same argument can also be used to prove that the same holds if the Hormander 
condition is satisfied only at a neighborhood of sup I or of any to G /. 

Example 16. Let M — (—2,2) x (— I, +00), with coordinates (x,y), and consider the vector 
fields 

((y + l)(l-x 2 ),-x) (x,h(x)(y + 2)) 

/( x ' v) = ,, =^yz 9{x, y) - 



yf{y + l) 2 {l-x 2 ) 2 +x 2 ' v / x 2 + h(x) 2 (y + 2) 2 
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^S5m=Rfo^i s 



iff 



Figure 1. The two vector fields of Example H6l with h(x) — ce 1 -*' 2 for x G [— 1, f] 



where h : [—2, 2] — > K is a smooth cutoff function such that supp/i C [— 1, 1], h > and /i(0) = f 
(see Figure [TJ. Fix < e < 1, C > 16 and let (/>, ip : [0, 1] — > M. be two smooth functions such 
that 



<Kt) 



H0<t<e, 
if 2e < t < 1, 



lW) 




if < i < 2e, 
if 3e < t < 1, 



and such that <f> is nonincreasing while ip is nondecreasing. Finally, consider the time-dependent 
system on M specified by the vector fields /*(x, y) = <f>(t)f(x,y) and g t (x,y) = ip(t)g(x,y), 
t € [0, 1]. We will show that {/*, g f } satisfies the Hormander condition for t G [0, e], but that the 
value function associated with the family {f t ,g t }te[o,i] 1S n °t lower semicontinuous. 

We start by showing that f(p) and g(p) are transversal for any p = (x, y) G M , implying the 
Hormander condition for {/*,<7*}, £ G [0, e]. If x G (—2,-1] U [1,2), then, by definition of h, 
g(p) = (1,0) is clearly transversal to /(p). On the other hand, if x G (—1,1) \ {0} and g(p) is 
parallel to /(p), a simple computation shows that h(x) < 0, which is a contradiction. Finally, for 
x = 0, it is clear that g(p) = (0, y + 2) and /(p) = (y + 1, 0) are never parallel. We remark that 
this implies also that the value function p £ , induced by controls defined on [0,e], is a distance 
equivalent to the Euclidean one. In particular, |pi — pi\ < 2p £ (pi,p2) for any pi,p2 G M. 

Fix now q' = (1, 0). The set of points from which q' is reachable using only / is exactly O q / = 
{(1,3/): y > -1}- Let then q G (-1,0) X {0} be such that p £ (q , (-1, 0)) < imin peC1<j , p e (q ,p). 
In order to show that p\{qo, •) is not lower semicontinuous at q', consider any sequence {g n }nGN C 
(1/2, 1) x {0} such that q n — > q' . By continuity of p e and the fact that — q„ — > (— 1, 0), we can 
always assume that, up to subsequences, p e (qo, ~Qn) < 5 niin pe ci , p e {qo,p)- 

Since y* = for t > 2e, if u G L 1 ([0, 1],R 2 ) is a control steering the system from qo to g', the 
control m| [o,2e] steers the system from qo to some p G O q i. Exploiting the fact that pi e > p £ by 
monotonicity of ip, this implies that 



(15) 



pi(qo,q') > min 

pec „' 



Pe(qo,P) > 2p e (go,-5n 
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Let now u G L 1 ([0,l],R m ) be the control constructed as follows. From time to e, w|[o, e ] is 
the minimizer of p E steering the system from q to — q n . Then, u|( e3e ) = and, after this, the 
control acts only on /' for time t G [3e, 1], steering the system from — q n to q n . Hence, 



(16) \q n -(-q n )\ = 



el 



el 



= C / \u(t)\dt. 

J3e 



/ u(t)f(x(t),y(t))dt 

he 

Since \q n - (-q n )\ < 2, C > 16/\q - q'\ > &/p e (qo,q'), and by (JTSJ, it holds that 

f 1 13 

pi(?0)?n) < / \u(t)\dt = p e (qo,-q n ) + ^\q n - (-<?«)! < ^Pi(<?o,<?')- 

Taking the liminf as n — > oo shows that pi(9oj •) is not l.s.c. at q' . 

3.3. Estimates on reachable sets. In this section, we concentrate on a particular class of 
time-dependent systems. Namely, let {/i,...,/ m } be a bracket-generating family of smooth 
vector fields, /o be a smooth vector field, and consider the time-dependent system 



(17) q = Yuift ft = ( e ~ t/o )*fi(q), q G M, u = («i 



n- 



Here, (e - *-^ )* is the push-forward operator associated with the flow of /( 

As we will see in the next section, this class of systems arises naturally when dealing with 
control systems that are afHne with respect to the control. Observe, in particular, that from the 
bracket-generating property of {/i, . . . , f m } it follows immediately that the time-dependent fam- 
ily {(e~*^°)*/i, . . . , (e -t ^°)»/ m } satisfies the strong Hormander condition, as per Definition [TT] 

Before proceeding with the estimates of the reachable sets, we need to define a suitable approx- 
imation of system (|17p . Namely, fix a system of privileged coordinates (in the sub-Riemannian 
sense) at q for {/i,.-.,/ m }- Assume that fo(q) / 0, and let s G {l,...,r} be such that 
ord q (/o) = — s. In this case, there exist, in coordinates, an homogeneous vector field /g^, of 
weighted degree — s, and a vector field /q S , of weighted degree > — s + 1, such that 

(18) z*fo = .fo S + tf- s - 
In particular, it holds that /q~ s ^ near z(q) = 0. 

Remark 17. The fact that ord 9 (/o) = — s is not equivalent, in general, to /o G A s near q. In 
particular, if the growth vector is non-constant around q, from ord ? (/o) > — s it does not follow 
that /o G A s . For example, consider the sub-Riemannian control system on R 2 with (privileged) 
coordinates (x,y), defined by the vector fields d x and xd y , called the Grushin plane. Outside 
{x = 0}, the non-holonomic degree of these vector fields is 1, while, on {x = 0}, we need one 
bracket to generate the y direction, and thus it is 2. Hence, if y ^ the vector field yd y is never 
in A near (0,y), but ord( ,jj) (yd y ) = ord (0iS) (y) + ord( ,y)(9 y ) = 0. 

However, due to the properties (0 and the fact that A s is a module, the converse is always 
true. Namely, if /o G A s , then ord g (/o) > — s. 

For any smooth vector field /, let (ad 1 f )f = [/ ,/] and {&d e f )f = [/o, (ad^Vo)/], for any 
£ GN. We recall (see for example [17]) that the following Taylor expansion holds 

(19) ( e - t/o )*/-E^( ad '/o)/- 

Since ord q ( fj) > — 1, by (|5]) we have that ord g ((ad fo)fj) > —£s — 1. Then, using the decomposi- 
tion (fT5|) , for any i > 0, there exists, in coordinates, an homogeneous vector field Ff of weighted 
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degree — is, and a remainder r e of order > — is — 1, such that 
(20) z*[(ad e f Q )f 3 ]=F' + r e . 

Definition 18. The homogeneous series approximation at q of ft, 1 < j < to, associated with 
the privileged coordinates z, is the vector field with coordinate representation 

(2D g=E* 






where g = L r_1 / ;5 J an d J" is the non-holonomic degree of {/i, . . . , / TO } at q. The approximated 
time-dependent control system is then defined as 

rn 

(ATD) 8 = X>(*)^G*)- 

If a system, in some system of privileged coordinates, coincides with its homogeneous series 
approximation , we will say that it is series homogeneous. 

The homogeneous series approximation encodes the idea that the time t is of weight s = 
— ord q (fo)- This is a consequence of the fact that, due to the expansion (fH?|) . t allows to 
build brackets of /o with the fjS. In this sense, the homogeneous series approximation is a 
generalization of the nilpotent approximation. 

We are now ready to state the main theorem of this section. 

Theorem 19. Let {/i, . . . , f m } satisfy the Hormander condition, let z = (zi, . . . , z n ) be a system 
of privileged coordinates at q £ M for {/i, . . . , fm}- Then there exist C, T, £o > such that, for 
any e < £o and an V <?' £ Ti-T(q,£), setting s — — ord g (/o) it holds 

(22) \ Zi (q')\<c(s Wi +eT^ if Wi <s, 

/ \wt-l 

(23) \zi{q')\<Ce[e + T-*) ifw t >s. 



Moreover, if the system is series homogeneous, then it holds the stronger estimate 

(24) \zi{q')\<Ce Wi if w t < s. 

To prove this theorem we need the following proposition, estimating the difference between 
(HU and (TXTDl . 

Proposition 20. Let {/i, . . . , / m } satisfy the Hormander condition, and let z — (z\, . . . , z„) be 
a system of privileged coordinates at q £ M for {/i, ...,/,„}. For T > and it £ L 1 ([0,T];R m ) 7 
let 7(-) anrf <y(>) &e the trajectories associated with u in (|17[) and (|ATD[) . respectively, and such 
that 7(0) = 7(0) = a. T/ien i/iere exist C,e$,Tq > 0, independent of u, such that, if t < To and 
J a \u\ ds = e < £q, and setting s = —otdq(fo) it holds 

(25) \z i { 1 {t))-z i m))\<Ce(e + r*) Wi , i = l,...,n. 

Remark 21. This proposition generalizes Proposition [B] In fact, in the sub-Riemannian case, 
since /o = 0, if, for t > 0, the curve 7 is associated to u £ L 1 ([0,t],R m ), it is associated also to 
u T (-) = ju(j-), for any r > 0. Thus, since /J" |u T | ds = f Q \u\ ds = e, ((25|) reduces to 

|«i( 7 (t)) - Zi(?Kt))\ < limC(e'"-+ 1 + r£ r ) = Ce w * +1 . 
Finally, assuming that u satisfies the hypotheses of Proposition [51 i.e., that |u| = 1, we get t = e. 
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Proof. Let z(j(-)) — x(-), z(-y(-)) = y(-), and ||z|| = YH=i \zi\ x ^ Wi - We mimic the proof of 
Proposition 7.29 in [8]. The first step is to prove that there exists a constant C > such that 
||x(t)||, ||y(i)|| < Ce for t and £ — J \u\ds small enough. We prove this for ||x(t)||, the same 
argument works also for ||y(i)||. 

In z coordinates, the equation of the control system (ITT]) is, 

m 

Xi(t) =22uj(t)(zi)*fj(~f(t)), i = l,...,n. 
3=1 
Due to the fact that z*f!j = z*fj + 0{t), uniformly in a neighborhood of q that ord 9 (zi) = Wi 
and that ord g (/j) > — 1, we have that there exist To and C > such that \(zi)*f]j(q)\ < 
%\(zi)*fj(q)\ < CHx^ir*- 1 , for any t < T . Thus we get 



(26) li.wifCC^K-wiiMOir- 1 . 

As in the proof for the sub-Riemannian case, choosing N sufficiently large, so that all N/wt are 
even integers, and setting |||z||| = (X)"=i \zt\ N ' Wi )x we get a norm equivalent to ||z||. Deriving 
with respect to time and using (j2"6")l we get ^|||x(t)||| < C^ =1 |uj(i)|. Finally, by integration, 
equivalence of the norms, and the fact that x(0) = z(q) = 0, we conclude that ||x(£)|| < Ce. 

Now we move to proving (|2"5|) . By construction of ([ATDp and the Taylor expansion of /* , for 
any I < g = [ r ~ 1 /s\ , there exist homogeneous polynomials h- of order wi —£s — l and remainders 
Tj i of order larger than or equal to Wi — is, such that we can write 

(^)*/i = E^+4)+^ e+1 )> 

£=0 l - 

(**)4 = E^4 

1=0 

Here, the O is intended as 1 1 and is uniform in a compact neighborhood of the origin. Then, 
m / q 

Xi(t)-in(t) = Y^ u i (*) ( 1_ , I 

j=l ^ £=0 



(t)-fc(t) = ^ tti (*) ^ ^ (^W - 4(2/) + 4(x)) + 0(^ +1 ) 



where Q l j ik are homogeneous polynomial in x and y, of order Wi — w k — is — 1. We observe that, 
if ^ — Wfe — is — 1 < 0, then Qj ik = 0. Thus, for sufficiently small ||a;|| and \\y\\, we have 

|^ ifc (z,y)| < C(\\x\\^- w "- is -^ + + \\y\\^- Wk - es - 1)+ ), \r%{x)\ < C\\x\\^-^ + . 

Here we let (£) + = max{£, 0}, for any (£l. Using the inequalities of the first step, taking t < T 
sufficiently small, and enlarging the constant C, we get 

\ii(t)-yi(t)\ 

<C\u{t)\(J2j( Y. \xk(t)-y k (t)\e Wi - Wk - es - l + e (wi - is)+ )+te +1 
^ e=o ' \w k <wi-is ' 

<C\u(t)\(j2 te ( E \x h (t)~y h (t)\s w '- Wh - 1 +s^-^ + )+t s+1 \ 
\ e =o \w h <wi ' ' 
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In the last inequality we applied the change of variable Wk >->■ Wh — £s in each of the sums. 

We can integrate the system by induction, since it is in triangular form. For Wj = 1, since 
[wi — £s) + — for any £ > 1, the inequality reduces to 

\xi(t) - yi(t)\ < C\u{t)\ [JVeC™*-^ +*e+ 1 j < C|u(t)|(e Wi +t). 
Wo / 

Here we enlarged the constant C. Thus, integrating the previous inequality, we get \xi(t)— yi(t)\ < 
Ce(e m +t) <Ce(e + t*) m . 

Let, then, Wi > 1 and assume that \xn{t) — Uh(t)\ < Ce(e + t~) Wh for any Wh < u>i. To 
complete the proof it suffices to show that \xi(t) — Vi(t)\ < C\u(t)\(e + t^) Wi , since (l2~5l) will 
follow, as above, by integration. Thus, we have, enlarging again the constant C and taking t 
sufficiently small, 

I*i(t)-W(*)l 

<C\u(t)\fj2t e ( J2 (e + tiy h e w *~ Wh +e^- es A+te +1 
(27) V e=o \w h <wi / 

< C\u(t)\(j2t e ( J2 t^e^- w -+s^- es A+t^ 
^ e=o \w h <wi J 

If t < e s , from (|2"T)) it is clear that \xi(t) — Vi{t)\ < C\u(t)\e Wi . Here we used the fact that 
q + 1 > Wi/ s. On the other hand, if e < t 1 ^ , it holds 



Q 

=0 \w h <w 



t t (t) - m\ < cnt)\ (j2 ( J2 t^ +e+ ^ K +t e+} ^ 11 +^ +i )<c\u(t)\t^ 



Putting together these two estimates, we get that \xi(t) — yi(i)\ < C\u(t)\(e Wi +t =' ) < C|u(i)|(£+ 
t~) Wi , completing the proof of the proposition. □ 

Proof of Theorem] 191 We start by claiming that (|M|) implies (j2"2")l . In fact, if 7 : q ~» q' is the 
trajectory associated in <jTTJ) to a control u G L 1 ([0,T],R m ), and 7 is the trajectory associated 
with the same control in the homogeneous series approximation (|ATD[) . with 7(0) = q, it holds 

ki(?')l<Ni(7(T))| + |^(7(T))-^(7(T))|. 

Thus, by Proposition I2U1 the claim is proved. 

Hence, from now on we assume our system to be in the form (jATDp . Let us define, for 
1 < j < n and < a < r, the vector helds tp? as 



T-F £ 
£=0 



^ -^ «'J ; 



where i*l £ are defined in ([20)1 . We do not explicitly denote the dependence on time, to lighten 
the notation. Observe that, if a = g, then, by (|2lT) . ip" = /■. 

We claim that, letting a;( Q )(-) be the trajectory associated with a control u G L 1 ([0, T],R m ) 
in system (|TD|) with {ip" , . . . , ip^} as vector fields, then, for some constant C > and any 
i G {1, . . . , n} and a > 1, it holds 

(28) \zF><X)-zf- 1 \T)\<l° , ix, ifWl " aS ' 
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In fact, due to the homogeneity of the F- , proceeding as in the proof of Proposition [301 we get 
that for Wi < as it holds 



\x^(t)-x^(t)\<c\u(t)\j2t £ E i4 a) (*)-4 a " 1) (*)ie"'- ,0h - 1 - 



By induction on 1 < w, < as, this proves the first part of the claim. On the other hand, if 
Wi > as, it holds 



\if\t) - ±t~ l \t)\ < C\u(t)\( J2 1" E l4 0) (*) " x { r l \t)\e Wi - Wh - X + re-— 1 

\ £=0 w h <Wi 

Then, again by induction over w„ we get that \x[ a) {T) ~ x[ a ~ l) (T)\ < CT a e w ^ as . Finally, the 
claim follows considering the cases T < e s and T > e s . 

Due to the fact that rf = fj, by Theorem it holds \x\ 0) (T)\ < Ce w *. Thus, applying $Zj$ 
and enlarging the constant C , we get 



, ( ,'. = Or )i 5:|..:-m-^. m | + |^ m |< ^ + rl ^^ 



This proves (|2"5jl and (|24p , completing the proof of the theorem. □ 

We end this section by showing that the estimate (|2"5)l is sharp, at least in some directions. 
Indeed, for a system which is series homogeneous at q in some privileged coordinates z, and 
satisfies the hypotheses of Theorem [T51 it holds that z*((ad fo)fj) is an homogeneous vector 
field of weighted degree — sk — 1. Thus, since et < e(e + t~) sk , the following proposition shows 
that (|23p is sharp in this direction. The proof is an adaption of an argument from [11] . 

Proposition 22. Let {/l, •••,/m} satisfy the Hormander condition. Let, moreover q G M, 
i £ {1, . . . ,m} and k > 0. Then, for any coordinate system y at q, there exist T, e$ > such 
that, for any e < eg and t < T there exists a (|TD|I - admissible curve 7 : [0, t] — > M, with 0(7) < e, 
and such that 

l/( 7 (t)) = et k dy((ad k f )f 3 (q)) + 0(et k+1 ) as et -+ 0. 

Proof. Let t, n > be fixed, and define u £ L 1 ([0,T],R m ) as Ui(r) = 77, Uj(t) = for j 7^ i, 
r S [0,i]. Then, fix any $ e C fc ([0, 1]) such that $«(0) = $«(1) = 0, for < i < k. Thus, by 
integrating by parts and the fact that ^(e ~ t ^°)*g — (e~*^ )*(ad(/o)5 r ), we get 

* (fc) (V0(e- T/o ),/i(g)dr = t k J^ir/tKe-^)* ((ad fe /o)/i) (g) dr, 

for any £ andq. This implies that the flows generated by & k \ T /t){e~ T f°)*fi a,ndt k $( T /t)(e~ T f°)* ((ad fc /o)/i 
coincide. Using the series expansions of the chronological exponential and (e^*^ ),, see [5J Section 
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2.4], there holds, then, 



p-t rn „t 

&$ V$ (fc) (V0%-M(e- T/o )*/;dT = e^/ t?*<*> (V0(e" rA )./i *■ 
Jo 7=1 Jo 

r^ fc *(V0(e- rA )*((ad fc /o)/i) dr 



= e3ql5 / 
Jo 

= c^| 7 ? t fc+1 $( S )(e-'^)* ((ad fc / )/i) ^ 

= c5^y ?/ t fc+1 *(,) ((ad fc /o)/i + 0(*)) rfs 
/•*+i/ , „,ifc* \ /. _i_ /mwfc+2 



= H + T^ 1 (adVo)/i + Ofar") 

Finally, considering any coordinate system and letting e = r/t, this completes the proof. □ 

4. CONTROL-AFFINE SYSTEMS 

In this section we apply the results of Section [3] to control- affine systems. Let {fi, ■ ■ ■ , f m } 
be a bracket-generating family of vector fields, /o be a smooth vector field, called the drift, and 
consider the control-affmc system 

m 

(D) g = /<>(?)+ £) «i/«(g), ? eM > »=K-.«m)ei m , 

»=i 

The assumption on {/i, . . . , / m } to be bracket-generating, is called strong Hormander condition 
for (0. 

An absolutely continuous curve 7 : [0,T] — > Af is (fDj-admissible if 7(f) = /o(7(^))+/u(t)(7(*)) 
for some control u € L 1 ([0,T],R m ). Its cost is defined as 

Cfoil) = min||u|| L i ([0 , T]iRm) , 

where the minimum is taken over all controls u such that 7 is associated with u. Then, we define 
the two value functions we are interested in as 

pT(q,q') = inf{c( 7 ) : 7 : [0,T"] -» M is ©-admissible, 7 : q ~* g' ', T" < T}, 

p*°(q, q') = inf{c(7): 7 (|D|) -admissible and 7 : g ^ q'}. 

It is clear that p^ilil') \ P {lil') as T — >■ +00, for any g, g' G M. Moreover, we observe 
that, pj!{q, e'^°g) = for any < t < T. Finally, the reachable sets with respect to these value 
functions, from any q G M and for s, T > 0, are 

1Z f T ° (q, e) = {q 1 G M : p f T ° (g, q') < e}, 7e /o (g, e) = {g' G M : /° (g, g') < e}. 

4.1. Connection with time-dependent systems. Applying the variations formula (see [5]), 
system (iDJ) can be written as a composition of a time-dependent system in the form (|17p and of 
a translation along the drift. Namely, for any u G L 1 ([0,T],R m ), it holds 

pT / "i v pT rn 

(29) e^ / f/o+53 «<(*) /<) di = eT/ ° ° ^ / 53 U 'W ( e ~ t/o )*/* *• 

J o V i=1 / Jo i=1 

We call time- dependent system associated with (|Dl) the following control system, 

771 

(30) q = Y. u ^ e ~ tf °)*^)' « eM ' u=(u 1 ,...,« m )eE m 
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Observe that, since diffeomorphisms preserve linear independence, the strong Hormander con- 
dition for (|D]), implies that {(e~^°)*/i, . . . , (e~ t: ' )*fm}te[o,+oo) satisfies the strong Hormander 
condition, defined in Definition 1111 

Exploiting these facts, we can prove the following. 

Proposition 23. Assume that (|L7|) satisfies the strong Hormander condition. Then, for any 
T > 0, the functions p^,p*° : M x M — >• [0, +oo) are continuous. Moreover, letting c?sr be the 
sub-Riemannian distance induced by {/i, . . . , f m }, for any q, q' £ M it holds 

p f T °(q,q') < mm T d S R(e tfo q,q'), p fo (q,q') < mmd SR (e tfo q,q'). 

Proof. The continuity of the two functions, and the fact that p^(q,q'), p^°(q,q) < dsR(q,q), for 
any q,q £ M, follows from the same arguments used in Theorem IT21 adapting Lemmata [TBI and 
ITU to the system (|D|). In particular, one has to consider (T, J 7 , £) \-t e'*° o E-r ?{£) instead of 

To prove the second part of the statement, we let, for any t £ [0, T), 

(p t (p) = inf{c(7) : 7 : [t, T'\ -> M is ©-admissible, 7 : p ~* g', T' < T}. 

It is clear that, as above, it holds iptip) < dsB.(p,q')- Moreover, we observe that p^{q,e t ^ q) = 
for any < t < T, and hence that for any such t it holds 

f#(q,q') < Vt{e tfo q) < ds R (e t/o g, «')• 
Taking the minimum for < t < T, proves the inequality regarding pjl. To complete the proof 
it suffices to observe that p^°(q, q') < Px(q, q') for any T > 0. □ 

We point out that in system ()D|) . as in time-dependent systems, the existence of minimizers 
is not assured. Moreover, this lack of minimizers is possible even if they exist for the associated 
time-dependent system, as the following example points out. 

Example 24. Consider the following vector fields on R 3 , with coordinates (x, y, z), 

fi(z,y,z) = d x> f 2 (x, y, z) =d y + xd z . 

Since [/1, f 2 ] = d z , {/1, f 2 } is a bracket-generating family of vector fields. The sub-Riemannian 
control system associated to {fi,f 2 } on R 3 corresponds to the Heisenberg group. 

Let now f = d z be the drift. Since [fi,d z ] = [f 2 ,d z ] = it holds that (e~*^°)»/i = f\ and 
(e _ *-' )*/2 = f 2 . Hence, the associated time-dependent system is actually not time-dependent. 
Thus, by (gg), a curve 7: [0,T] -> R 3 is ([SR|-admissible for {/i,/ 2 } if and only if 7 (-) = e /o o 7 (-) 
is (|D])- admissible. As a consequence of this, for any q £ R 3 and any e > 0, 

n f °(q,e) = \Je tfo oB SK (q,e). 

t>Q 

As pointed out in Section [2J minimizers for the sub-Riemannian system exist between any 
pair of points in 2?sr(<7, e), if £ is sufficiently small. Let us show that, for any point in TZ^°{q,e) 
with positive z coordinate, we have an explicit minimizer, while for the others there exists no 
minimizer. Without loss of generality we can consider q = 0. Then, since e* *°(x', y', z') — 
(x' t y', z' + 1'), every point (x, y, z) £ 1Z^ (0, e) with z > 0, can be reached optimally considering 
the sub-Riemannian minimizing curve between the origin and (x, y, 0) rescaled on time z. 

If, instead, z < 0, we cannot construct any sub-Riemannian trajectory from to (x, y, z — t), 
t > 0, with cost < dsR(0, (x,y, z)). This is due to the fact that the minimizing trajectories in 
Heisenberg group are the lifts of arcs on the plane (x, y), spanning area equal to the z coordinates, 
and that \z — t\ = — z + t > \z\. Since, by Proposition [531 /° /o (0, (x, y, z)) < dsR(0, (x 7 y,z)), this 
implies that there exists no minimizer for p'°(q, (x,y, z)). 
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4.2. Estimates on reachable sets. In this section we apply Theorem \T§[ in order to obtain 
results in the spirit of Theorem [7] and Corollary [5J First, we need the following definition. 

Definition 25. The point q is said to be regular with respect to the drift /o, if q' i-4 ordq<(/o) is 
locally constant at q. 

The main result of this section are the following local regularity estimates for pf° . We cannot 
expect anything global, since in general the sets TZ^°(q,e) are noncompact. 

Theorem 26. Assume that (iDJ) satisfies the strong Hormander condition, and let q be regular 
with respect to the drift /o- Assume, moreover, that z = (zi,...,z n ) is a system of privileged 
coordinates at q for {/i, . . . , f m }, such that Z*fo = d Zk , for some 1 < k < n. Then, there exist 
Tq,Eq,C > such that 

1 dist (ztf), z{e^ T ^q)) < p f T ° (q, q') < C dist (z{q'), z{e^ T ^q)) ^ , q' G K f T ° (q, s), 

where, for any x G M™ and A C W l , dist(x, A) — inf^g^ \x — y\ is the Euclidean distance between 
them and r is the degree of non-holonomy of {/i, . . . , f m } at q. 

Let us define the following sets, for parameters rj > and T > 0. We remark that Box (77) is 
defined as in (flU)) and that {d Zi }f =1 is the canonical basis in W 1 . 

Srfa)= (J (W+Box^A 

U T (v) = Box (77) U (J {z G R" : z k = £, \z t \ < n w ' + 77^ for w t < s, 1 ± k, 

0<i<T 

and \zi\ < 77(77 + ^) Wi ' 1 for w { > s}, 
Ut(v) = Box (77) U (J {z G M" : z k = £, \z t \ < n w ' for w, < s, 1 ± k, 

0<£<T 

and \zi\ < 77(77 + ^) Wi ~ l for Wj > s}. 

As in the sub-Riemannian case, Theorem [2H1 is a direct consequence of some estimates on the 
shape of the accessible sets, contained in the following. 

Theorem 27. Assume that (iDJ) satisfies the strong Hormander condition, and let q G M be 
regular with respect to the drift /o. Assume, moreover, that z — [z\, . . . , z n ) is a system of 
privileged coordinates at q for {/1, ...,/,„}, such that z*fo = d Zk , for some 1 < k < n. Then, 
there exist C, £o,T > such that 

(31) i^c^McW, /o r£ < E0 „,T<T„. 

Here, with abuse of notation, we denoted by TU£ (q, e) the coordinate representation of the reach- 
able set. In particular, 

Box f^ej n {z k < 0} C K%(q,s) D {z k < 0} C Box(Ce) n {z k < 0}. 

Moreover, if the system is nilpotent, it holds 

(32) ^'.jc.jMcta M £ < e „«,r<r . 

In order to prove Theorem 1271 we need the following lemma. 
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Lemma 28. Let z = (z\, . . . , z n ) be a system of privileged coordinates at q G M. Then there 
exist C,Eo,Tq > such that, for any q' G TV£ (q, £o) for £ < £o and T < To, and such that 

(i) for any t < s , ord g /( t ) / = -s, where q'{t) = e~ tfo {q'), 

(ii) dzk(fo(z(q'))) 7^ 0, for some k with Wk = —s, 

it holds that, if u G L 1 ([0, T], K m ) is a control steering the system ([D)l from q to q' , with \\u\\i = e, 
then 

T < C(e s +max{z fc ( g '),0}). 

Proof. For any n > 0, let 7 be the trajectory associated with u G L 1 ([0,T],R m ) in the system 
(|Dj) . and satisfying 7 : q ~~» q'. Let 7 be the trajectory associated with u and starting from q, in 
the time-dependent system (j3"0]) . Thus 7(f) = e*^° 07(f) and p(q,j(T)) < e. 

Recall that, for any vector field g and point p £ M, it holds that z k (e Tg (p)) — z k (p) = 
/ dzk(g(e t9 (p))). Thus, by the mean value theorem, there exists r G [0,T] such that 

(33) z k (q') = z k ( 7 (T)) = rdz fc (/o(e T/o (7(T)))) + %(7(r)). 

Since e T ^°('j(T)) — e~( T ~ T ^°{q'), by hypothesis (ii) and the smoothness of /o, there exist 
T ,Ci > 0, independent of 7, such that dz k (f Q (e T f"(^(T)))) > d for T < T . Hence, by 
Theorem \T§\ (since w k = s), there exist C2, £ > such that, if e < e and T < To, 

|z fc ( 7 (T))|+max{z fc ( g '),0} C 2 (e a + Tg) +max{z fc (g / ),0} 

Ci d ■ 

Since the constants are independent of 7, taking C = C2/C1, £0 < min{T ,£, (C — 1)/C 2 } 
completes the proof. □ 

Proof of Theorem \2T\ The first inclusion in pip follows from Theorem [7] and Theorem 1231 In 
fact, combining them, we have that, for any e < Eq and any T > 0, 

E T (^e)c |J B SR (e t f°q,e)cl4°(q,e). 

^ ' 0<t<T 

To prove the second inclusion, we let q' G TZ^(q,e). Fix any 77 > and consider a control 
u G L 1 ([0,T],R m ) such that its associated trajectory 7, in the system |D]), satisfies 7 : q ~-> 5' 
and c/ (7) < £ + ?y. We distinguish two cases. First we assume that Zk{q') < 0. In this case, by 
Lemma [23 it follows there exists C, £o,Tb > such that if T < To and e < £0, then T < Ce s . 
Moreover (l29l) implies that e~ T f°(q') G 1Zt (q,e). Then, enlarging the constant C, Theorem H9l 
yields 

\z t (q')\ = \ Zl (e' Th {q'))\ < C(e m +eT^) < Ce w \ if w, < s and i ^ k, 

\z k (q')\ <T +\z k (e- T ^(q'))\ <T + C(e + T*) S <Ce s , 

\ Zl (q')\ = |^(e- T ^(Q'))| <Ce(e + T 1 ^- 1 <Ce w % if Wi > s. 

Here, we used the fact that, for any p G M, from z*/o = d Zk , it holds Zi(p) = Zi(e~ T ^°(p)) and 
\dz k (fo(p))\ = I- Thus, g' C Box(C£) cn(Ce). 

On the other hand, if Zk(q') > 0, Lemma 1251 yields that T < C(£ s + Zfc(g'))- Then, applying 
again Theorem Q30 we get 

\zi(q')\ < C(e w ' +ez k (q')^), iiwi<s and i ^ k, 

\z k (q')\ <T + Ce s , 

\z t {q')\ < Ce(e + z k {q')i) m ~\ if w % > s. 



This proves that 9'cII (Ce), completing the proof of ^1 
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To prove Q32I) it suffices to use the same argument as above, applying the result on nilpotent 
systems in Theorem ItUl □ 

Remark 29. Theorem |2"T1 suggests that the behavior of the system (|D]l. when moving in the 
direction — /o, is essentially sub-Riemannian. However, although this is true locally in time, it is 
false in general. For example, consider the Euclidean plane endowed with a rotational drift, i.e., 
such that {e*^°(g)} te (o i +oo) is diffeomorphic to S 1 for any q ^ 0. Then, p'° (q, e - *-* (q)) = for 
any t > and thus we can move in the direction — /q for free. 



Proof of Theorem \26l Since every norm on R™ is equivalent, dist(z(q'), [0, T]d Zk ) is equivalent to 

a(q')= V \Zi(q')\+ min \z k (q r ) - 1\. 

l<2<n 

Thus, to complete the proof it suffices to prove that it holds C^ x a{q') < p^(q,q') < Ca(q'Y^ r . 

By Theorem [TH S T (C _1 e) C K^(q,e) C n T (Ce) for any e < e . The first inclusion is 
equivalent to the fact that, for every e < £o such that Ca(q') < e r ', one has p^(q, q') < s. From 
this follows that p^(q,q') < C 1 / T a{q') 1 / r . The same reasoning applied to the other inclusion 
proves that 

\zi(q')\ < C(p f T °(q,q') Wi + P f T °(Q,Q')T^) if w t < s, i ± k, 
mhi\z k tf)-t\<Ct#(q t <{)', 

\z t {q')\ < C(p f T °(q,q>) w * + p f T °(q,q')T^) if w , > s. 

Clearly, this implies that a(q') < Cp^'(q,q'), for some larger constant, completing the proof of 
the theorem. □ 
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